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Chapter  1 
Introduction 


This  is  the  final  report  on  Contract  No.  F33615-88-C-1792  for  Wright-Patterson  Air 
Force  Base  over  the  period  August  1988  through  May  1991.  The  overall  goals  and 
accomplishments  are  summarized.  Some  of  the  highlights  of  the  major  technical 
findings  are  presented  in  the  various  chapters.  More  detailed  information  can  be 
found  in  the  reports  and  papers  listed  below. 

1.1  Statement  of  Problem 

The  overall  goal  of  this  program  is  to  assess  the  use  of  the  polarization  response  of  a 
scattering  object  as  added  information  along  with  modern  SAR  techniques  in  order 
to  identify  or  discriminate  a  radar  target.  In  particular,  the  program  has  three  major 
tasks.  First,  to  provide  high  quality  calibrated  measurements  on  a  set  of  basic  targets 
for  the  complete  scattering  matrix  and  both  backscatter  and  bistatic  scattering  using 
The  Ohio  State  University  ElectroScience  Laboratory  (OSU-ESL)  compact  range. 
Second,  to  measure  specific  scale  model  targets  in  the  OSU-ESL  compact  range  in  the 
frequency  range  2-18  GHz  and  full-scale  targets  at  the  Big  Ear  range  at  frequencies 
below  2  GHz.  In  both  cases,  coherent  magnitude  and  phase  signal  versus  target 
orientation,  frequency  and  polarization  are  collected.  Third,  to  analytically  model 
the  basic  targets  in  order  to  provide  insight  into  the  characteristics  of  the  individual 
scattering  mechanisms.  The  modeling  effort  is  based  on  the  Uniform  Geometrical 
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Theory  of  Diffraction  (UTD).  Computer  codes  for  the  various  hasic  targets  have 
resulted. 

The  basic  targets  of  concern  are  flat  plates,  dihedrals,  trihedrals,  top-hats,  “cake 
pans,”  bruderhedrals,  bi-cones,  and  toruses.  Some  of  these  objects  are  chosen  for 
their  usefulness  as  calibration  targets.  It  is  also  assumed  that  real  targets  of  interest 
as  they  appear  in  Synthetic  Aperture  Radar  (SAR)  are  composed  of  similar  objects. 

1.2  Accomplishments 

Modern  UTD  techniques  have  been  used  to  study  basic  shapes  applicable  to  target 
identification.  First  order  terms  related  to  simpler  physical  optics  approximations, 
and  higher  order  terms  have  been  developed  to  model  the  structures.  The  higher 
order  terms  are  important  to  obtain  more  complete  lower  frequency  behavior.  In 
some  cases  correct  amplitude  information  in  lower  cross  section  regions  is  totally  lost 
without  this  type  of  component.  Potential  characteristic  “arcs”  in  the  SAR  images 
also  can  be  associated  with  these  terms. 

Analysis  of  the  rectangular  and  triangular  dihedral  has  been  completed.  A  com¬ 
puter  code  that  provides  patterns  and  swept  frequency  information  has  been  written 
and  validated  against  method  of  moments  and  the  above  mentioned  measurements. 
It  can  be  used  to  validate  measurements,  identify  physical  mechanisms  in  the  re¬ 
sults,  perform  parametric  studies  in  designing  dihedrals  of  various  edge  shapes,  and 
conduct  component  studies  of  more  general  teal  targets. 

The  analytic  techniques  have  been  extended  to  top-hats  and  the  concave  shape 
of  a  “cake  pan,”  which  can  represent  a  wheel  hub.  Bruderhedrals  and  bi-cones  have 
also  been  studied.  Preliminary  results  for  torus  shapes,  representing  tires,  have  been 
obtained.  Computer  codes  have  been  written  to  validate  the  analysis.  In  addition, 
simpler  formulae  have  been  obtained  to  help  provide  insight  into  possible  target 
identification  characteristics. 
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Simple  flat  plate  stnictures  have  been  extensively  studied  to  obtain  amplitude 
and  frequency  behavior  over  wide  bands  and  angles.  A  much  more  complete  de¬ 
piction  for  scattering  from  the  corners  has  been  obtained.  In  addition,  the  higher 
terms  that  propagate  across  the  plates  from  corner  to  edges  and  corners  to  corners 
have  been  studied.  The  results  have  been  validated  against  other  analyses  and  high 
resolution  measurements. 

The  analytic  results  obtained  have  been  used  to  obtain  a  better  understanding 
of  SAR  images.  They  have  also  been  used  to  help  design  calibration  targets  in  SAR 
scenes.  The  analysis  and  measurements,  help  clarify  phenomenology  pertinent  to 
target  identification. 

A  radar  measurement  system  for  obtaining  Ultra- Wideband  (UWB)  scattering 
data  over  the  47-  to  2000-MHz  band  has  been  developed  and  used  on  full  size  ground 
vehicle  targets  up  to  70  ft  long.  Data  have  been  gathered  on  several  targets  over  a 
wide  range  of  radar  look  angles. 

A  technique  for  calibration  of  UWP  ranges  using  dihedrals  and  trihedrals  has 
been  developed.  The  technique  has  been  used  for  comparison  measurements  of  three 
diagnostic  scattering  systems. 

A  UWB  linear  imaging  process  has  been  developed  and  implemented.  Diagnostic 
images  have  been  produced  for  ground  vehicles  (both  models  and  full  size),  as  well 
as  for  several  of  the  fundamental  scattering  targets  which  have  been  analyzed  by  this 
study.  It  has  been  shown  that  the  images  have  sufficient  resolution  and  sensitivity 
to  resolve  higher  order  scattering  mechanisms,  which  are  the  subject  of  our  advances 
in  GTD  analysis. 

An  image  windowing  algorithm  which  permits  the  extraction  of  the  pattern  and 
spectrum  of  a  single  scattering  mechanism  has  been  developed  and  demonstrated. 
The  characteristics  of  several  higher  order  scattering  mechanisms  have  been  verified 
^sing  this  algorithm. 
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Details  of  the  above  accomplishments  and  more  are  snmma  ’  ,ed  in  tliis  report 
and  in  the  following  reports  and  papers. 

The  following  reports  have  been  generated  relative  to  this  contract: 

•  “Far  Zoi.e  Scattering  from  Rectangular  and  Triangular  Dihedrals  in  the  Prin¬ 
cipal  Plane,”  by  N.  Akhter  and  R.  J.  Marhefka,  Technical  Report  721198-1, 
February  1989.  (Ako  a  M.Sc.  Thesis) 

•  “Ultra-Wideband  Ra  'ar  Imaging  and  the  Diagnosis  of  Scattering  Centers,” 
by  S.  M.  Scarborough,  Technical  Report  721198-2,  September  1990.  (Also  a 
M.Sc.  Thesis) 

•  “An  Assessment  of  Three  Experimental  Radar  Scattering  Measurement  Sys¬ 
tems  Operating  Below  2  GHz,”  by  J.  S.  Gwynne,  Technical  Report  721198-3, 
March  1991.  (Also  a  M.Sc  Thesis) 

•  “High  Frequency  Electromagnetic  Scattering  by  Polygonal  Structures,”  by 
L.  P.  Ivrissimtzis  and  R.  J.  Marhefka,  Technical  Report  721198-4,  March  1991. 
(Also  a  Ph.D.  Dissertation) 

•  “'^urrent  Near  the  Vertex  of  a  Perfectly  Conducting  Angular  Sector,”  by  T.  J. 
Brinkley  and  R.  J.  Marhefka,  Technical  Report  721198-5,  June  1991.  (Also  a 
Ph.D.  Dissertation) 

The  foil  ving  on!  papers  have  been  presented  in  association  with  this  contract: 

•  “Polarimetric  ISAR  Imaging  to  Identify  Basic  Scattering  Mechanisms  of  Ge¬ 
ometrically  Complicated  Objects,”  G.  Dural,  J.D.  Young  and  D.L.  Moffatt, 
1989  IEEE  APS-URSI  Meeting,  San  Jose,  Cabfornia. 

•  “Comparison  of  Three  Low-Frequency  Measurement  Techniques  and  Calcula¬ 
tions  for  Small  Dihedral  Backscattering,”  J.S.  Gwynne  and  J.D.  Young,  1989 
IEEE  APS-URSI  Meeting,  San  Jose,  Cabfornia. 
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•  “Image  Processing  for  Diagnosis  of  Dispersive  Scattering  Centers,”  S.M.  Scar¬ 
borough  and  J.D.  Young,  1989  IEEE  APS-URSI  Meeting,  San  Jose,  California. 

•  “A  Study  of  Scattering  from  Edge  Waves  Using  High  Resolution  Radar  Im¬ 
ages,”  J.D.  Young,  1990  IEEE  APS-URSI  Meeting,  Dallas,  Texas. 

•  “Ultra-Wideband  Diagnostic  Target  Imaging,”  J.D.  Young  and  S.M.  Scarbor¬ 
ough,  First  Los  Alamos  Symposium  on  Ultra- Wideband  Radar,  May  1990. 

•  “UTD  Scattering  from  Dihedrals,”  by  R.  J.  Marhefka  and  N.  Akhter,”  5th 
Annual  Review,  Applied  Computational  Electromagnetics  Society,  Naval  Post¬ 
graduate  School,  Monterey,  California,  March  21-24,  1989. 

The  following  written  papers  have  been  submitted  in  association  with  this  con¬ 
tract: 

•  “Current  Near  the  Vertex  of  a  Perfectly  Conducting  Angular  Sector,”  by  T.  J. 
Brinkley  and  R.  J.  Marhefka  submitted  to  IEEE  Transactions  on  Antennas 
and  Propagation. 

•  “Far  Zone  Field  Scattered  by  a  Perfectly  Conducting  Angular  Sector,”  by  T.  J. 
Brinkley  and  R.  J.  Marhefka  submitted  to  IEEE  Transactions  on  Antennas  and 
Propagation. 

1.3  Overview 

The  following  report  discusses  both  the  measurements  and  analysis  performed  on 
this  contract.  Chapter  2  discusses  the  UTD  analysis  of  a  two-dimensional  dihedral. 
The  application  of  UTD  to  a  top-hat  and  bi-cone  is  given  in  Chapter  3.  Chapter  4 
presents  a  newly  developed  diffraction  coefficient  for  the  scattering  from  a  corner 
of  a  plate.  The  higher  order  terms  on  a  plate  are  briefly  discussed  in  Chapter  5. 
An  ultra-wideband  imaging  technique  for  characterizing  the  frequency  and  pattern 


behavior  of  local  scattering  renters  is  presented  in  Chapter  6.  The  cnlibrntion  and 
performance  of  three  measurement  systems  and  the  various  results  of  simple  targets 
obtained  are  outlined  in  Chapter  7.  Chapter  8  gives  the  conclusions. 
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Chapter  2 

Two-Dimensional  Dihedral 

2.1  Introduction 

UTD  is  used  to  analyze  backscattet  and  bistatic  scattering  Irom  an  arbitrary  angled 
dihedral.  The  solution  contains  interactions  of  up  to  third  order,  which  provides  ex¬ 
cellent  results  for  dihedrals  of  60°  or  greater  and  for  sizes  larger  than  one  wavelength 
in  extent.  Parallel  ray  formulation  of  the  UTD  is  used  here  for  the  far  zone  problem. 
This  means  that  only  edge  diffractions,  their  images,  and  various  combinations  of 
double  diffractions  are  needed.  This  is  a  very  efficient  and  accurate  means  of  ana¬ 
lyzing  the  far  zone  problem.  It  avoids  any  integrations  and  nonphysical  “imposed” 
edge  conditions. 

The  UTD  interactions  that  have  been  included  in  the  analysis  and  contained  in 
the  2-D  computer  code  are  shown  in  Figure  2.1.  The  single  interaction  fields  included 
are  simply  edge  diffraction  (D)  terms.  The  second  order  interaction  fields  are  the 
reflected  -  diffracted  (RD),  diffracted  -  reflected  (DR),  and  diffracted  -  diffracted 
(DD)  (double  diffracted)  terms.  The  third  order  interaction  fields  are  the  reflected 
-  reflected  -  diffracted  (RRD),  diffracted  -  reflected  -  reflected  (DRK  f,  reflected  - 
diffracted  -  reflected  (RDR),  reflected  -diffracted  -  diffracted  (RDD),  diffracted  - 
diffracted  -  reflected  (DDR),  and  diffracted  -  reflected  -  diffracted  (DRD)  terms. 
Triple  diffractions  have  not  been  included. 
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Figure  2.1:  The  interactions  considered. 


2.2  Basic  Dihedral  UTD  Mechanisms 


The  main  beam  of  the  pattern  is  formed  by  the  fields  from  RD  and  the  DR  UTD 
terms  along  with  the  single  diffracted  (D)  field  from  the  wedge  formed  at  the  junc¬ 
tion  between  the  two  intersecting  strips  providing  the  necessary  fields.  The  double 
diffraction  (DD)  across  the  face  of  the  dihedral  can  also  be  important. 

For  angles  derived  from  Figure  2.2  for  bistatic  scattering  from  a  general  angled 
dihedral,  the  diffracted  term  from  the  junction  is 


D\R  and  RDi  result  due  to  reflection  from  strip  2  and  the  terms  D^R  and  RDi 
result  due  to  reflection  from  strip  1.  The  RD  and  DR  fields  are  given  by 
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where,  the  reflection  coefficient,  R^  h  =  q:l. 

It  is  noted  that  the  term  D\R  is  present  in  the  region 
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Also,  RDi  and  RD2  are  present  in  the  regions  given  by  Equations  (2.4)  and  (2.5), 
respectively,  with  ^  replaced  by  <f)'. 

These  terms,  in  genera],  look  quite  complicated.  In  some  special  cases,  they  can 
be  reduced  to  simpler  formulas.  Consider  the  case  of  a  =  90°.  Using 
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Close  to  backscatter,  with  kl\  large,  the  second  term  is  dominant.  Thus  the  above 
expression  may  be  written  as 


rrO,ft+/iZJi  _  rri  ^k/4  ^  ^  ^ 

^dominani  —  ^ s,h^  ^  I  2  I 

I  sin  [2W,  sin  (^)  sin  (^)]  1 

[  2W,6in(^)sin(^)  J 


(2.8) 


For  the  backscatter  case  with  strips  of  equal  length,  from  Equations  (2.3)  to  (2.5) 
and  (2.8) 

0<*<V4-  (2-9) 

and 

’'/2-  (2-'«) 

Notice  that  for  <f>  =  45°  this  result  agrees  with  the  well  known  PO  result.  Also,  ob¬ 
serve  that  there  is  a  slope  discontinuity  in  the  dominant  UTD  (RD  and  DR  terms) 
at  ^  =  45° .  The  DD  field  across  the  face  of  the  dihedral  balances  this  slope  discon¬ 
tinuity. 

The  specular  from  the  individual  strip  faces  is  produced  by  D  from  the  ends  of 
the  strips  and  the  DR,  RD  and  DRD  terms  between  the  strips. 


2.3  More  Multiple  Plate  UTD  Terms 

Depending  on  the  angle  of  the  dihedral  and  on  the  bistatic  angle,  various  terms  can 
be  important.  The  RRD  and  the  DRR  terms  are  present  when  the  angle  of  the 
dihedral  is  less  than  90°.  Of  course,  these  two  terms  are  most  significant  when  the 
diffraction  is  along  a  reflection  shadow  boundary  RSB.  Although  the  RDR  term  is 
present  for  all  angles  of  the  dihedral,  its  significance  is  greatest  when  the  angle  of  the 
dihedral  is  less  than  90°  and  the  diffraction  is  along  a  reflection  shadow  boundary. 

For  the  backscatter  case  with  strips  of  equal  lengths,  based  on  simple  reflection 
concepts,  the  order  of  interactions  required  to  get  the  pattern  between  the  strips 


is  the  largest  integer  whose  magnitude  does  not  exreed  Thus,  for  n  —  60°, 

third  order  terms  are  required  to  get  the  pattern  between  the  strips.  Figure  2.3 
shows  how  only  the  third  order  terms  for  a  dihedral  of  60°  combine  to  give  the 
correct  result  between  the  strips.  Figure  2.4  compares  the  Method  of  Moments 
results  with  that  of  the  UTD  computer  code  (with  all  terms  present)  for  a  dihedral 
of  60°.  Notice  that  the  discontinuity  at  <^  =  60°  in  Figure  2.3  and  Figure  2.4  is  due 
to  the  triply  diffracted  field  which  has  not  been  included  in  the  computer  code.  Our 
formulation  indicates  that  for  a  =  45°,  fourth  order  interactions  are  required.  Since 
the  UTD  computer  code  does  not  take  into  account  the  fourth  order  interactions, 
we  expect  no  agreement  of  the  results  between  the  strips  in  this  case.  This  is  shown 
in  Figure  2.5.  However,  even  for  this  case,  notice  the  excellent  agreement  for  ^ 
greater  then  50°.  Thus,  for  the  terms  present  in  the  computer  code,  the  pattern 
would  degrade  gracefully  between  the  strips  for  a  less  then  60°.  As  the  angle  of  the 
dihedral  becomes  smaller,  more  terms  are  required  to  give  the  correct  result  betw»*en 
the  strips  and  a  “cavity”  effect  results. 

2.4  Results  of  Rectangular  Dihedral 

In  this  section,  the  theory  and  computer  code  are  validated  by  comparing  the  UTD 
results  to  2-D  method  of  moments  results.  In  order  to  calculate  the  backscatter 
results,  a  slight  bistatic  angle  is  used  to  avoid  numerical  limit  concerns  on  the  com¬ 
puter.  Method  of  Moments  is  chosen  to  validate  the  results  as  an  independent  check 
since  both  are  two-dimensional.  Fairly  large  sized  dihedrals  can  be  calculated  using 
method  of  moments  in  the  2-D  case.  In  3-D,  limits  on  computer  memory  restrict  the 
size  of  dihedral  that  can  be  analyzed  using  method  of  moments.  A  9"  (22.86  cm) 
2-D  dihedral  is  analyzed  at  2  GHz  and  10  GHz  meaning  the  dihedral  is  1.52  and  7.62 
wavelengths  on  a  side,  respectively.  All  the  patterns  are  normalized  to  dB  relative 
to  a  square  meter,  that  is,  a  two-dimensional  cross  section. 
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Figure  2.3:  Only  third  order  terms  give  a  continuous  result  for  9"  60°  rectangular 
dihedral  at  10  GHz  for  a  horizontally  polarized  field  (backscatter). 
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Figure  2.5;  Backscatter  field  from  a  9"  45°  rectangular  dihedral  at  10  GHz  for  a 
horizontally  polarized  field. 
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For  the  vertically  polarized  incident  field,  the  backscatter  from  a  90°  dihedial, 
Figure  2.6  shows  the  singly  diffracted  (D)  term,  DR,  and  RD  terms.  Figure  2.7 
shows  the  D,  RD  and  DR  terms  combined,  the  DRD  term,  and  the  DD  term.  This 
is  all  that  is  required  to  get  the  complete  result  for  the  90°  dihedral.  Figure  2.13 
shows  the  result  obtained  by  the  D,  DR,  RD,  DRD,  and  DD  terms  combined. 

From  here  on,  the  UTD  results  contain  all  of  the  terms  discussed  above.  For  the 
backscatter  case,  from  Geometrical  Optics,  a  main  beam  is  expected  when  the  angle 
of  the  dihedral  is  such  that  is  an  integer.  This  is  clearly  evident  in  Figures  2.4, 

2.5,  and  2.7.  Also,  as  the  angle  of  the  dihedral  diverts  from  a  which  makes  an 

integer,  the  main  beam  breaks  up.  Figures  2.8  and  2.9  show  backscatter  results  for 
a  =  85  and  95  degrees,  respectively,  with  the  incident  field  horizontally  polarized. 
As  expected,  the  main  beam  of  the  pattern  for  a  =  90°  breaks  up  as  a  is  decreased 
or  increased. 

The  backscattered  field  for  the  2-GHz  case  is  shown  in  Figure  2.10  for  horizontal 
polarization  and  in  Figure  2.11  for  vertical  polarization.  The  10-GHz  backscatter 
case  is  shown  in  Figures  2.12  and  2.13  for  the  horizontal  and  vertical  polarizations, 
respectively.  Patterns  for  a  fixed  bistatic  angle  of  90°  are  shown  next.  The  2-GHz 
case  is  shown  in  Figure  2.14  for  horizontal  polarization  and  in  Figure  2.15  for  vertical 
polarization.  The  10-GHz  patterns  are  shown  in  Figures  2.16  and  2.17  for  horizontal 
and  vertical  polarization,  respectively.  In  all  cases,  the  agreement  is  excellent.  The 
small  dihedrals  have  slightly  more  disagreement  because  the  higher  order  interactions 
are  more  important.  It  is  assumed  even  better  agreement  and/or  smaller  sized 
dihedrals  could  be  achieved  and/or  analyzed  by  including  more  interactions.  For 
engineering  accuracy,  however,  this  would  not  seem  necessary. 

2.5  Conclusion 

The  UTD  formulation  gives  excellent  results  for  the  rectangular  dihedral  in  the 
principal  plane.  For  backscatter,  the  number  of  interactions  required  to  get  the 
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Figure  2.8:  Backscatter  from  85“  rectangular  dihedral  at  10  GHz  for  horizontally 
polarized  held. 


Figure  2.9:  Backscatter  from  95“  rectangular  dihedral  at  10  GHz  for  horizontally 
polarized  field. 
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igure  2.10:  Backscatter  at  2  GHz  for  horizontally  polarized  field. 
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Figure  2.11:  Backscatter  at  2  GHz  for  vertically  polarized  field. 
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Figure  2.14;  90°  bistatic  field  at  2  GHz  for  horizontally  polarized  field. 


Figure  2.15:  90°  bistatic  field  at  2  GHz  for  vertically  polarized  field. 
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complete  pattern  for  a  dihedral  of  angle  a  is  found  to  be  the  integral  vabie  of 
For  bistatic  scattering,  the  number  of  interactions  required  wiU  be  less  than  or  equal 
to  that  for  the  backscatter  case.  Since  the  computer  code  developed  for  this  analysis 
contains  interactions  of  up  to  third  order,  our  results  are  found  to  be  valid  for 
backscatter  for  dihedrals  of  60°  or  greater.  This  solution  is  expected  to  degrade 
gracefully  between  the  strips  for  dihedrals  of  less  than  60°.  In  general,  as  the  angle 
of  the  dihedral  diverts  from  a  which  makes  the  an  integer,  the  main  beam 

amplitude  breaks  up  rapidly. 
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Chapter  3 

Top-Hats  and  Bi-Cones 

3.1  Introduction 

Due  to  the  localized  nature  of  high  frequency  scattering,  the  solution  to  complex 
structures  may  be  modeled  as  a  superposition  of  scattered  fields  from  simpler  shapes. 
The  object  of  this  research  has  been  to  find  far  zone  backscatter  and  bistatic  scattered 
fields  from  some  of  the  most  common  shapes.  This  -equires  second  and  higher  order 
interactions  t  .  get  accurate  results  over  a  wide  range  of  angles.  For  these  shapes, 
the  most  important  higher  order  interactions  have  been  determined.  Where  possible, 
the  expressions  will  be  simplified  to  show  the  connection  between  the  UTD  results 
and  the  more  classical  PO  approximations. 

The  two-dimensional  dihedral,  studied  in  detail  in  [1],  is  one  of  the  simplest 
of  such  shapes.  The  result  of  the  two-dimensional  dihedral  may  be  modified  by 
curvature  and  caustic  corrections  to  get  results  for  the  top-hat,  bi-cone,  bruderhedral, 
cake  pan,  et  cetera.  Other  shapes  that  have  been  studied  are  the  donut  and  trihedral. 

Ray  techniques,  like  the  Geometrical  Theory  of  Diffraction,  predict  infinite  fields 
at  the  caustics.  The  occurrence  of  these  caustics  requires  the  erroneous  field  infini¬ 
ties  be  “corrected”  with  caustic  correction  functions  which  blend  the  wide-angled 
predictions  smoothly  into  the  correct  value  at  the  caustic.  Caustic  corrections  of  the 
second  and  higher  order  are  required  to  get  accurate  results  for  the  shapes  mentioned. 
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Caustic  corrections  for  the  second  and  liigher  order  terms  have  been  investigated  in 
the  context  of  UTD. 

The  solutions  to  the  top-hat  and  the  bi-cone  are  discussed  in  Sections  3.2  and  3.3, 
respectively.  The  curvature  modifications  and  the  caustic  corrections  to  the  two- 
dimensional  solution  are  shown  to  give  appropriate  results  when  compared  with  the 
body  of  revolution  results,  even  for  very  electrically  small  targets. 

3.2  Top-Hat 

3.2.1  Introduction 

The  2-D  dihedral  result  of  Reference  [1]  may  be  used  to  get  3-D  results  for  a  wide 
variety  of  problems.  The  solution  of  the  triangular  dihedral  in  the  principal  plane 
is  also  shown  in  [1],  The  solution  to  the  top-hat  problem  is  discussed  in  this  report. 
A  few  modifications  are  required  to  the  UTD  2-D  analysis  of  the  dihedral.  Analytic 
results  are  validated  against  body  of  revolution  results.  Also,  for  backsca*ter,  a 
dominant  UTD  result  is  shown  to  be  valid  in  all  regions  and  is  reasonably  simple  to 
program.  The  parameters  for  the  top-hat  are  shown  in  Figure  3.1. 

3.2.2  Modifications  Required 

Basically  rim  caustic  correction  and  spread  factors  are  required  to  modify  the  2-D 
dihedral  solution  to  give  the  solution  for  the  top-hat.  In  Reference  [1],  the  most  dom¬ 
inant  terms  for  the  90°  dihedral  are  found  to  be  diffraction,  diffracted- reflected  (DR), 
reflected- diffracted  (RD),  and  diffracted-reflected-diffracted  (DRD).  The  dominant 
contribution  from  these  terms  will  be  used  to  get  the  result  for  the  top-hat. 

3.2.3  Theoretical  Background 

In  addition  to  the  background  information  present  in  Reference  [1],  we  need  in¬ 
formation  on  the  rim  caustic  correction.  For  the  backscatter  case,  when  the  two 
diffractions  are  occurring  from  the  same  point  (as  in  the  DRD  case),  the  double 
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Figure  3.1:  Top-hat  geometry. 

diffraction  formulation  can  be  simplified.  The  simplified  equations  are  showi.  below, 
in  Section  3.2.4. 

3.2.4  Rim  Caustic  Correction 

Details  of  the  single  diffraction  rim  caustic  correction  are  shown  in  References  [2] 
and  [3].  Basically,  a  half  rim  is  associated  with  each  stationary  point.  The  stationary 
points  being  the  two  diffraction  points  outside  the  caustic  region.  To  correct  for  the 
rim  caustic  the  regular  diffraction  coefficient  is  modified  by  a  transition  function. 
This  result  goes  to  the  diffraction  point  solution  outside  the  caustic  region,  making 
it  a  uniform  solution.  The  :  'ffied  diffraction  coefficient  is  given  as 

where 

0  =  0^=  (3.2) 


25 


r,(t/)  =  i/,(f/)y^e-'("-5)  (3.3) 

and 

/.(t/)  =  »[^-ff,(£/)  +  jV,(£^)].  (3.4) 

Ji  and  Hi  are  the  first  order  Bessel  and  Struve  functions,  respectively.  Also, 


U  =  karie  •  (f  +  r') . 


(3.5) 


The  transition  function  is  Te  {U).  k  is  the  wave  number,  a  and  fie  are  the  radius  of 
curvature  and  the  normal  to  the  edge  at  the  point  of  diffraction,  and  f  and  f'  are 
the  directions  to  the  source  and  the  receiver,  respectively.  The  transition  function 
has  small  and  large  argument  forms  given  as: 
for  \U\  «  1 


for  \U\»  I 
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It  is  important  to  note  when  U  is  negative,  the  negative  branch  of  the  square  root 
is  selected. 


3.2.5  Double  Diffraction 


Details  of  double  diffraction  formulation  used  here  are  presented  in  [4].  For  the 
backscatter  case  when  the  two  diffractions  occur  from  the  same  point  and  Op  =  a^, 
the  double  diffraction  equations  for  P„  and  Ph  can  be  simplified.  We  have 
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PA^,n,d)  = 
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for  a  TM  plane  wave  (with  the  electric  fiehl  parallel  to  the  wedge).  For  a  TE  plane 
wave  (with  the  electric  field  perpendicular  to  the  wedge), 


PH{^,n,d)  = 


(3.9) 


These  equations  are  used  to  get  the  dominant  contribution  from  the  DRD  terms  of 
the  side  and  top  of  the  top-hat.  All  the  remaining  equations  remain  unchanged  and 
are  given  in  Reference  (1). 


3.2.6  Dominant  Result 

As  mentioned  above  the  dominant  result  is  obtained  from  the  diffracted,  DR,  RD, 
and  DRD  terms.  We  will  use  the  exact  expressions  for  all  the  terms,  except  for  the 
case  of  DRD  where  only  the  discontinuous  term  is  considered. 


3.2.7  D,  RD,  DR,  and  DRD  Terms 


The  expressions  shown  below  are  for  backseat  ter  from  a  top- hat  of  parameters  given 
by  Figure  3.1.  These  expressions  are  valid  for  —90°  <  <f>  <  90°  unless  otherwise 
stated.  Of  course  by  symmetry,  the  result  for  the  other  side  may  be  obtained. 

The  diffracted  terms  from  the  disk  for  the  right  side  is  given  by 
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where  R^_h  =  and 


V..  =  £  ■  F.  field  perpendicular  to  the  pattern  plane 

*’  I  /T  ,  H  field  perpendicular  to  the  pattern  plane. 

The  diffracted  terms  from  th«.  ^op  of  the  cylinder  from  the  right  side  is 
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The  diffraction  reflection  (£>|  refers  to  diffraction  from  the  edge  of  the  disk  and 
R  refers  to  the  reflection  from  the  cylindrical  surface)  terms  are  given  by 


t/o  5  0  <0<  arctan  y 

Y2cos(^)  [l\ 

and  (f?2  refers  to  diffraction  from  the  top  and  R  refers  to  the  reflection  from  the 
disk) 
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Y  2  cos  (^) 

For  the  DRD  term,  only  the  dominant  term  is  retained  for  the  joined  wedge 
solution.  The  dominant  DRD  term  from  the  side  of  the  top-hat  is 
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b  +  (cos(^)  -  1)(/  -f  p’’) 
and  ($,n,<i)  is  given  by  Equation  (3.8)  and  (3.9).  The  dominant  DRD  term 
from  the  right  side  of  the  top  of  the  top-hat  is 
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The  dominant  DRD  term  from  the  left  side  of  the  top  of  the  top-hat  is 
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3.2.8  Peak  of  the  Top-Hat 

Along  the  axis  of  the  top-hat,  the  cotangent  terms  of  the  diffraction  coefficients 
that  tend  to  “blow  up”  combine  to  give  the  dominant  result.  Not  surprisingly,  the 
dominant  result  obtained  is  exactly  as  one  would  predict  from  physical  optics. 

As  expected,  the  reflection  boundary  cotangent  terms  of  the  diffraction  from  the 
upper  face  of  the  disk  (obtained  from  Equations  (3.10)  and  (3.11))  give  the  area  of 
the  disk.  The  reflection  boundary  cotangent  terms  from  the  diffraction  from  the  end 
cap  of  the  upper  face  of  the  cylinder  (obtained  from  Equations  (3.13)  and  (3.14))  give 
the  area  of  the  end  cap.  The  incident  boundary  cotangent  terms  of  the  curved  face 
from  DR,  RD  and  from  DRD  (obtained  from  Equations  (3.16),  (3.21)  and  (3.22)) 
give  the  negative  area  of  the  end  cap.  Thus  the  apparent  “infinities”  combine  to 
give  a  result  equal  to  the  area  of  the  disk  minus  the  area  of  the  end  cap  of  the 
cylinder  plus  the  area  of  the  end  cap  of  the  cylinder  with  the  phase  accounting  for 
twice  the  electrical  distance  of  the  cylinder.  The  RCS  along  the  axis  of  the  top-hat 
is  appropriately  given  as 

=  47r 

3.2.9  Results 

Figures  3.2  and  3.3  show  the  results  obtained  for  the  horizontal  and  vertical  polar¬ 
ization,  respectively,  by  the  UTD  computer  code  and  the  Body  of  Revolution  code 
at  10  GHz  for  a  =  2  in.,  b  =  4  in.,  and  h  =  4  in.  top-hat.  The  results  compare  very 
well. 


7r(6^  —  a^)  7ra 


-t- 


-jk2h 


(3.23) 


3.3  The  Bi-Cone 

3.3.1  Introduction 

A  solution  for  the  bi-cone  is  sought  by  using  the  UTD  2-D  solution  of  the  dihedral 
discussed  in  Reference  [1].  The  parameters  for  the  bi-cone  are  shown  in  Figure  3.4. 
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Figure  3.4:  Bi-Cone  geometry. 

3.3.2  Modifications  Required 

Basically,  caustic  correction  and  spread  factor  modifications  are  required  to  the  2-D 
dihedral  solution  to  get  a  result  for  the  bi-cone.  Two  types  of  caustics  need  to  be 
corrected.  One  is  the  rim  caustic  which  has  been  used  in  the  solution  to  the  top 
hat  in  Section  3.2.4  with  details  given  in  [2]  and  [3].  The  other  is  the  curved  surface 
caustic.  Since  diffraction  is  a  local  phenomena,  it  cannot  be  expected  to  account  for 
a  caustic  of  a  global  phenomena.  However,  the  curved  surface  caustic  is  applicable 
if  the  two  ends  of  the  curved  surface  have  the  same  radii  of  curvature  (same  opread 
factors  for  the  two  diffraction  points).  The  diffraction  coefficient  from  each  of  the 
edges  “blows  up”  in  the  specular  direction  but  combines  to  give  the  correct  result. 
With  the  two  radii  of  curvatures  different,  the  diffraction  coefficients  from  the  two 
edges  have  different  spread  factors  and  the  diffraction  coefficient  cannot  combine  to 
give  the  correct  result. 
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3.3.3  Theoretical  Background 

In  addition  to  the  background  information  presented  in  Reference  [1]  and  in  Sec¬ 
tion  3.2,  we  need  information  on  the  curved  surface  caustic  correction. 

3.3.4  Curved  Surface  Caustic  Correction 

First  Order 


Details  of  the  single  diffraction  curved  surface  caustic  correction  are  presented  in  [2]. 
Basically,  the  term  of  the  diffraction  coefficient  which  “blows  up”  in  the  specular 
direction  is  corrected  by  the  inclusion  of  a  transition  function.  This  result  goes  to  the 
diffraction  point  solution  outside  the  caustic  region,  making  it  a  uniform  solution. 
The  specular  terms  from  the  two  diffraction  points  with  the  correction  are  given  as 


%(+jXa,) 


(3.24) 


The  angles  <f>  and  <!>'  are  measured  from  the  curved  surface,  Oe  is  the  radius  of 
curvature  of  the  diffracting  edge,  and 

X  =  if  cot  Qc  Jtan^ttc  (cos^  V -f  sin^  -f  ij  ^ 

•  sin/3o2cos  j  cos  j  .  (3.26) 


The  transition  function  is  given  as 


T(z)  =  T,(z)  =  -zM(l,2.5,-z) 


(3.26) 


(3.27) 


where 


T(2)  =  T2(z)=  -2ff(l,2.5,-z) 


C/(l,2.5,-z)=  -?M(l,2.5,-z)+^(-z)-ie- 


(3.28) 


(3.29) 
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Pliysirally,  T|  may  be  viewed  as  a  result  due  to  a  finite  rone  frtistum,  and  T2  due  to 
a  semi-infinite  cone  frustum.  Mathematically,  both  give  the  same  result  when  both 
edges  of  the  reflecting  surface  are  considered.  The  small  and  the  large  argument  of 


the  transition  function  are  the  following: 

for  \z\  <<  1, 

,  .,(22  22^  2"\ 

(3.30) 

(3.31) 

for  |2:|  >>  1, 

(3.32) 

(3.33) 

3.3.5  Higher  Order  Terms 

To  maintain  the  physical  optics  result  for  the  cone  frustum  in  the  specular  region, 
the  argument  of  the  transition  function  for  the  curved  surface  must  be  the  same  for 
DR,  RD,  and  DRD  terms.  However,  where  is  the  radius  of  curvature  at 

the  reflection  point.  The  formulation  remains  the  same  as  for  the  first  order  curved 
caustic  correction  discussed  above,  in  3.3.4. 

3.3.6  Dominant  Result 

As  expected,  th**  bi-cone  has  the  same  type  of  dominant  mechanisms  as  the  2-D 
dihedral  and  the  top-hat.  Thus,  the  dominant  result  is  obtained  from  diffracted, 
DR,  RD,  and  DRD  terms  for  the  bi-cone. 

3.3.7  D,  RD,  DR,  and  DRD  terms 

The  expressions  shown  below  are  for  backscatter  from  a  bi-cone  with  parameters 
given  in  Figure  3.4.  These  expressions  are  valid  for  —90°  <  <p  <  90°  unless  otherwise 
stated.  Of  course,  by  symmetry,  the  result  for  the  other  side  may  be  obtained. 
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The  diffraction  coefficient  from  the  top  on  the  right  side  is  given  hy 
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The  diffraction  coefficient  from  the  bottom  on  the  right  side  is  given  by 
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The  diffraction  reflection  terms  (D|  refers  to  the  diffraction  from  the  top  right 
edge  and  R  refers  to  the  reflection  from  the  curved  side  of  the  bottom  cone)  are 
given  as 


Ijlh  H  ^  l/ROi  ^  ij. 
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where 


T,  =  Te  ^jfcre  cos  (^  -  <^))  (3.39) 

and  {D2  refers  to  the  diffraction  from  the  bottom  right  edge  and  R  refers  to  the 
reflection  from  the  curved  side  of  the  top  cone) 
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A'  =  6  + 


cos 

cos  (I  -  *)  cos  (1) 


(3.46) 


(3.47) 


The  curved  surface  caustic  correction  for  the  double  diffraction  ensures  conti¬ 
nuity  of  the  fields.  For  the  diffracted  -  reflected  and  reflected  -  diffracted  fields, 
the  caustic  correction  is  required  on  the  term  associated  with  the  incident  shadow 
boundary.  Using  the  separate  wedges  formulation,  assuming  angles  for  the  double 
diffraction  are  measured  from  the  curved  surface,  the  terms  P»,h  (^p=2,  ^11=2)  and 
P»,h  (^p=3)^q=3)  insure  continuity  of  the  incident  shadow  boundary  of  the  DR  and 
RD  terms.  Thus,  the  transition  function  required  for  the  curved  surface  caustic 
correction  is  included  only  with  these  two  terms.  However,  using  the  joined  wedges 
formulation  and  assuming  that  angles  are  measured  from  the  end  cap,  the  term 
P»,h  (^p=t,^9=i)  insures  continuity  of  the  DR  and  RD  terms. 

For  the  DRD  term,  only  the  dominant  term  is  retained  for  the  joined  wedges 
solution.  The  dominant  DRD  term  at  the  top  of  the  bi-cone  (especially  around 
<t>  =  45°)  is 
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and  from  the  bottom  (especially  around  (j)  =  —45°) 
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where 
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3.3.8  Results 


(3.50) 
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Figures  3.5  and  3.6  show  the  results  obtained  for  the  horizontal  and  vertical  polar¬ 
ization  by  the  UTD  computer  code  and  the  Body  of  Revolution  code  at  10  GHz  for 
a  =  2  in.,  b  =  0.5  in.,  and  I  =  2\/L5  in.  bi-cone,  where  I  is  the  slant  length  of  one 
of  the  cones.  It  is  expected  that  by  including  the  diffracted- diffracted  field  the  UTD 
results  should  become  even  better. 


3.4  Conclusions 

It  has  been  shown  that  excellent  results  for  electrically  small  sizes  of  a  top-hat  and 

bi-cone  are  obtained  by  making  the  necessary  curvature  modifications  and  caustic 

corrections  to  the  2-D  dihedral  solution.  It  is  expected  that  with  similar  modifica- 

» 

tions,  solutions  to  the  bruderhedral  and  the  cake-pan  may  be  obtained.  The  solution 
to  the  donut  also  requires  caustic  corrections  of  first  and  second  order. 
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Figure  3.5:  Backscatter  field  from  a  bi-cone  for  horizontally  polarized  field.  The 
solid  line  is  UTD  solution  and  the  dashed  line  is  the  BOR  solution. 


^  (deg.) 


Figure  3.6:  Backscatter  field  from  a  bi-cone  for  vertically  polarized  field.  The  solid 
line  is  UTD  solution  and  the  dashed  line  is  the  BOR  solution. 
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Chapter  4 

Feir  Zone  Corner  Diffraction 


4.1  Introduction 

Many  local  characteristic  scattering  centers  of  ground  targets  can  be  modeled  by 
flat  plate  structures.  In  the  far  zone,  the  first  order  effects  show  up  in  images  as 
emanating  from  the  corners  of  the  plates.  The  UTD  treats  this  as  the  incident  field 
multiplied  by  a  diffraction  coefficient.  Comparisons  of  previous  corner  diffraction  so¬ 
lutions  with  measurement  based  on  ultrawideband  images  have  shown  discrepancies 
in  some  regions  of  space. 

The  report  “Current  Near  the  Vertex  of  a  Perfectly  Conducting  Angular  Sector,” 
by  T.  J.  Brinkley  and  R.  J.  Marhefka  explores  a  new  component  missing  in  previous 
results.  An  approximate  expression  for  the  current  on  an  angular  sector  is  derived 
from  the  exact  solution  [5].  It  is  used  to  find  a  new  far  zone  expression  for  the 
field  diffraction  by  a  vertex.  This  new  vertex  component  is  added  to  parts  based 
on  PO  and  fringe  currents  [6]  and  higher  order  local  edge  interaction  parts  [7].  The 
resulting  diffraction  coefficient  is  summarized  in  the  next  section. 

The  new  expression  for  the  corner  diffracted  field  is  verified  by  comparison  with 
measurements  where  the  contribution  from  a  single  vertex  has  been  singled  out  using 
time  domain  gating  techniques.  The  results  are  also  applied  to  large  flat  plates 
and  compared  with  excellent  agreement  to  measured  patterns.  The  advantage  for 
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target  identification  is  that  the  local  frequency  and  angular  behavior  of  vertex  type 
structures  can  now  be  more  accurately  predicted. 


4.2  Corner  Diffracted  Field  Expression 


The  field  diffracted  by  the  corner  of  an  angular  sector  is  written  as  the  sum  of  the 
contribution  to  the  scattered  field  from  the  different  currents  on  the  angular  sector. 
The  far  zone  field  scattered  by  an  angular  sector  is 
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(4.1) 
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(4.11) 
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The  various  angles  are  shown  in  Figure  4.1.  The  details  of  the  derivation  may  be 
found  in  [8],  [9],  or  [10]. 

The  contribution  from  the  edge  diffracted  current  of  edge  m  (m  =  1,2)  termi¬ 
nated  by  edge  m  (m  =  1,2)  is 
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The  functions  C*p,  C*p,  C^,,  5',  and  cos”*  have  been  defined  previously  in  Equa¬ 
tions  (4.6),  (4.7),  (4.8),  (4.12),  and  (4.10),  and  the  angles  are  defined  in  Figure  4.1. 
The  details  of  the  derivation  may  be  found  in  [9]  or  [10]. 

The  contribution  from  the  vertex  diffracted  current  is 
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incident  ray 


projection  of  the 
incident  ray 
onto  the  yg2*c 


projection  of  the 


Figure  4.1:  Edge  fixed  coordinates  on  an  angular  sector. 
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The  angles  are  defined  in  Figure  4.2.  Expressions  for  the  currents  used  in  Equa¬ 
tion  (4.19)  are  given  in  [11].  The  number  of  eigenvalues  and  the  number  of  terms 
used  in  calculating  the  eigenfunctions  are  =  7  and  N^j  =  9  to  calculate  J  ''‘’(f’) 
for  the  surface  integral  in  Equation  (4.18),  Nm  =  15  and  N^f  =  15  to  calculate 
J  )  for  the  line  integral  in  Equation  (4.18),  and  Ne„  =  25  and  iV,./  =  15  to  cal- 
culate  J  '‘’(t^)  for  the  line  integral  in  Equation  (4.18).  The  values  of  the  parameters 
used  in  all  of  the  calculations  are  r',  =  0.2A,  rj  =  0.3A,  =  0.5,  and  =  1.5 

which  have  been  found  to  give  fairly  accurate  results,  while  limiting  the  computa¬ 
tional  time  required.  The  integrations  in  Equation  (4.18)  must  be  done  numerically. 
A  simple  Riemann  sum  is  used  to  do  the  integrations  using  4  steps  for  the  radial 
part  of  the  surface  integral,  20  steps  for  the  angular  part  of  the  surface  integral,  and 
20  steps  for  the  line  integration. 
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Figure  4.2:  Radiation  integral  geometry  for  a  far  zone  field  point. 


The  contribution  from  the  edge  diffracted  current  of  edge  m  (m 
nated  by  the  other  edge  (m  =  2, 1)  is 
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The  details  of  the  derivation  may  be  found  in  (12,  7]. 

The  contribution  from  the  double  diffracted  current  terminated  by 
diffracting  edge  is 
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The  angles  are  defined  in  Figure  4.1,  and  the  details  of  the  derivation  may  be  found 
in  (7,  12]. 

4.3  Results 


To  illustrate  the  accuracy  of  the  new  corner  diffracted  field,  it  is  compared  with 
other  methods.  The  results  labeled  the  “New  Corner”  in  Figure  4.6  are  calculated 
using  the  contributions  from  the  physical  optics  current,  edge  diffracted  currents, 
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z 


edge  1 


Figure  4.3:  Vertical  triangle. 

the  double  diffracted  current,  and  the  vertex  diffracted  currents  as  given  in  the 
section  above.  The  parameters  used  in  calculating  the  contribution  from  the  vertex 
diffracted  current  are  =  15  and  N^j  =  15  to  find  J 7Ve„  =  25  and  N^j  =  15 
to  find  J  r',  =  0.2A,  =  0.3A,  20  equally  spaced  rectangular  approximations 

for  the  numerical  evaluation  of  the  linear  integral,  Nev  =  7  and  Ne/  =  9  to  calculate 
J  for  the  surface  integral,  and  4  equally  spaced  rectangular  approximations  in  the 
r'  part  of  the  surface  integral,  and  20  equally  spaced  rectangular  approximations  in 
the  part  of  the  surface  integral.  The  results  are  compared  with  measured  results 
and  the  previous  corner  diffraction  coefficient  based  on  the  contributions  from  the 
physical  optics  current  [8,  9, 10],  edge  diffracted  currents  [9, 10],  and  double  diffracted 
currents  [7,  12]. 

Measurement  data  are  processed  using  the  following  procedure  to  provide  a  ref¬ 
erence  solution  for  the  field  diffracted  by  the  corner  of  an  angular  sector.  Frequency 
swept  measurements  over  a  frequency  range  of  2  to  18  GHz  are  made  on  the  triangle 
shown  in  Figure  4.3  for  both  the  case  a  =  32",  a„  =  30°  and  the  case  a  =  32", 
Ov  =  45°  at  5°  degree  intervals  in  the  9  =  90°  pattern.  The  frequency  swept  data  at 
each  pattern  point  are  transformed  to  the  time  domain  to  obtain  the  band  limited 
impulse  response.  The  approximate  impulse  response  is  convolved  with  a  35-dB 
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Figure  4.4:  Impulse  response  polarized  co-pol  response)  of  the  triangle  in  Fig¬ 
ure  4.3  with  a=32"  and  a„  =  30°  at  =  90°  and  4>  =  145°. 

Kaiser-Bessel  smoothing  waveform  in  order  to  eliminate  Gibb’s  phenomenon  ring¬ 
ing  at  the  system  cut-off  frequency  of  18  GHz.  The  time  response,  typical  of  those 
obtained  for  pattern  points  in  the  region  95°  <  <f>  <  175°,  for  the  pattern  angle 
<!>  %  145°  is  shown  in  Figure  4.4.  The  two  first  order  mechanisms  are  labeled  El, 
which  corresponds  to  the  return  from  the  front  edge,  and  C3  which  corresponds  to 
the  return  from  corner  3.  The  higher  order  mechanisms  labeled  CE  and  CC  corre¬ 
spond  to  corner-edge  interactions  and  corner-corner  interactions.  The  path  length 
to  corner  3  is  sufRciently  different  from  the  path  lengths  of  the  other  mechanisms 
that  the  field  diffracted  by  corner  3  may  be  gated  out  in  the  time  domain  over  most 
of  the  pattern.  The  response  from  corner  3  is  isolated  in  the  time  domain  and  then 
transformed  back  to  the  frequency  domain  to  obtain  the  contribution  from  this  cor¬ 
ner  over  the  entire  frequency  band.  Further  details  of  the  time  gating  process  may 
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Figure  4.5:  Frequency  response  polarized  co-pol  response)  of  corner  3  of  the 
vertical  triangle  with  a=32"  and  =  30°  at  ff  =  90°  and  ^  =  145°. 

be  found  in  [13]  and  [14].  The  contribution  from  corner  3  as  a  function  of  frequency 
is  shown  in  Figure  4.5. 

The  measurements  may  then  be  used  to  determine  the  scattered  field  from  an 
angular  sector  as  a  function  of  the  pattern  angle.  From  the  frequency  dependence  of 
all  of  the  components  of  the  corner  diffracted  field  in  the  section  above  and  Figure  4.5 
(which  gives  the  cross  section  in  dBsm),  the  corner  diffracted  field  has  a  frequency 
response  that  varies  as  one  over  the  frequency.  The  corner  diffracted  field  is  therefore 
constant  for  a  given  look  angle  when  it  is  expressed  in  terms  of  dBSW.  The  roll-off 
at  the  upper  and  lower  limits  of  the  frequency  response  in  Figure  4.5  is  a  result 
of  the  time  gating  process,  so  only  the  data  in  the  center  of  the  frequency  range 
are  used.  The  measurements  are  averaged  (the  contribution  from  the  corner  should 
be  a  constant  when  it  is  expressed  in  dBSW)  over  50  evenly  spaced  samples  in  the 
region  of  6  to  14  GHz  to  obtain  the  measured  value  of  the  corner  diffracted  field. 
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Repenting  this  process  for  all  of  the  points  in  the  pattern  where  the  individual  corner 
contribution  may  be  gated  out  in  the  time  domain  results  in  the  measured  curves  in 
Figure  4.6.  The  results  show  that  the  additional  correction  to  the  corner  diffracted 
field  improves  the  high  frequency  calculations  considerably  for  this  pattern. 
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Figure  4.6:  Calculated  and  measured  contribution  from 
of  the  vertical  triangle  for  a  0  =  90®  pattern. 


Chapter  5 

Higher  Order  Plate  Terms 

5.1  Goals 

The  scattering  from  polygonal  plate  structures  is  composed  of  its  first  order  effects 
as  discussed  in  the  Chapter  4  as  well  as  higher  order  terms  that  travel  across  the 
plates.  These  higher  order  terms  show  up  as  arcs  in  the  SAR  images.  The  goal  of 
this  effort  is  to  understand  the  phenomenology  of  these  higher  order  mechanisms  in 
order  to  determine  their  frequency  response  and  pattern  behavior.  This  information 
is  intended  to  provide  insight  into  whether  these  mechanisms  are  robust  enough  to  act 
as  characteristics  identifying  tags  of  targets.  In  addition,  preliminary  determination 
of  diffraction  coefficients  for  double  diffraction  between  noncoplanar  skewed  edges 
has  been  accomplished.  These  types  of  diffraction  coefficients  are  needed  for  the 
evaluation  of  the  complete  patterns  of  trihedrals  and  dihedrals  in  three  dimensions. 
The  detuls  are  given  in  the  report  “High  Frequency  Electromagnetic  Scattering  by 
Conducting  Polyhedral  Structures,”  by  L.  P.  Ivrissimtzis  and  R.  J.  Marhefka.  Some 
relevant  results  to  this  contract’s  goals  are  outlined  here. 

5.2  Method 

The  phenomenon  of  higher  order  interactions  between  edges  in  flat  plate  structures 
for  monostatic  and  bistatic  cross  section  is  investigated.  Specifically,  a  ray  approxi- 
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mation  of  the  total  field  is  developed  starting  with  a  postulate  of  tlie  current  on  the 
scatterer.  The  surface  current  includes  a  primary  nonuniform  component  generated 
by  the  edges  of  the  structure  in  addition  to  the  Physical  Optics  constituent.  A  sec¬ 
ondary  nonuniform  current  excited  by  doubly  diffracted  fields  is  also  superimposed. 
The  asymptotic  reduction  of  the  related  surface  or  bne  integral  results  in  individ¬ 
ual  mechanisms  that  are  formulated  in  terms  of  corner-to-edge  and  edge-to-corner 
diffracted  rays. 

Edge  wave  interactions  are  examined  as  a  field  component  in  a  high  frequency 
formulation  of  the  scattering  by  edged  structures.  The  edge  wave  corner  launching 
diffraction  is  derived  via  reciprocity  from  the  approximation  of  the  diffraction  of 
dipole  excited  edge  waves.  Rediffraction  of  tip  launched  edge  waves  is  formulated  by 
replacing  the  tip  that  excites  them  with  a  fictitious  dipole  source  radiating  close  to 
one  of  the  edges  of  a  trihedral.  Pattern  and  time  domain  calculations  reveal  that  the 
estimation  of  the  edge  wave  interaction  is  adequate,  as  long  as  this  type  of  return 
does  not  coincide  in  time  with  corner  diffraction  or  doubly  diffracted  rays. 

The  canonical  problem  of  near  zone  double  diffraction  by  a  pair  of  coplanar 
skewed  edges  is  studied  by  employing  an  Extended  Spectral  Theory  of  Diffraction. 
A  closed  form  double  diffraction  coefficient  resulting  from  the  asymptotic  analysis 
is  derived,  and  may  be  employed  in  pattern  calculations  of  sources  radiating  in  the 
presence  of  conducting  flat  plate  objects.  The  solution  is  extended  to  the  case  of 
noncoplanar  skewed  edges  via  a  heuristic  modification  of  the  pertinent  transition 
function,  so  that  continuity  across  the  shadow  boundaries  of  the  edge  diffracted  rays 
is  warranted. 

5.3  Results 

The  various  types  of  the  higher  order  mechanisms  that  can  exist  on  a  plate  are 
illustrated  in  Figure  5.1.  It  depicts  a  24"  square  plate  tilted  at  a  45°  angle  with 
respect  to  the  horizontal  plane.  The  measured  images  are  presented  in  Chapter  6  in 
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Figures  6.5  and  6.6.  The  aspect  angle  is  scanned  from  —5°  to  95°  in  step  of  0.5°  wliile 
for  each  angle  the  frequency  is  scanned  from  2  to  18  GHz  in  steps  of  20  MHz,  The 
arcs  appearing  in  the  images  of  the  plate  correspond  to  higher  order  mechanisms. 

The  corner  diffracted  pulses  are  denoted  by  C1-C4  while  possible  edge  wave 
mechanisms  are  indicated  by  E1-E4.  The  letter  D  is  used  to  indicate  doubly 
diffracted  rays.  In  order  to  compare  the  calculated  and  measured  results  more  clearly, 
down  range  cuts  are  compared  in  the  complete  report  discussed  in  this  chapter.  Fig¬ 
ures  5.2,  5.3,  and  5.4  compare  measurements  to  calculations  for  the  6  =  30°  vertical, 
horizontal,  and  cross-polarized  components,  respectively.  Note  that  in  the  measured 
data,  the  two  returns  (Si,  S2)  are  from  the  mount  supports.  Unknown  mechanisms 
are  marked  with  question  marks. 

In  most  cases  good  agreement  is  obtained.  This  is  particularly  true  for  the  leading 
corners;  trailing  corners  disagreed  in  some  cases.  This  should  be  improved  with  the 
new  vertex  component  discussed  in  Chapter  4.  It  was  not  available  at  the  time  of 
this  report. 

To  observe  the  frequency  behavior  of  isolated  mechanisms,  several  of  the  sepa¬ 
rable  ones  have  been  time  gated  and  transformed  back  into  the  frequency  domain. 
Figures  5.5  and  5.6  show  the  B  =  30°  cases  for  edge  wave  mechanisms  marked  El 
and  E2,  respectively.  Excellent  agreement  between  calculations  and  measurements 
are  obtained  for  many  of  the  cases.  The  largest  discrepancy  is  associated  with  the 
double  diffraction,  D,  terms  across  the  plate.  To  illustrate  the  behavior  of  these 
terms,  the  time  range  between  2-3  nsec  has  been  gated  and  transformed.  These 
comparisons  are  shown  in  Figures  5.7  and  5.8.  Apparently,  the  frequency  character¬ 
istic  of  a  double  diffraction  mechanism  is  comparable  with  that  of  an  edge  wave.  As 
a  matter  of  fact,  in  the  high  frequency  domain,  theory  predicts  both  should  decay 
with  the  inverse  of  The  shape  of  these  measured  and  computed  curves  in  the 
frequency  domain  appear  very  similar. 
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The  discrepancies  in  the  curves  are  most  likely  attributable  to  the  fact  that 
these  higher  order  mechanisms  have  been  analyzed  based  on  the  simpler  corner 
diffraction  models.  It  will  be  interesting  to  see  if  the  more  advanced  model  described 
in  Chapter  4  can  be  used  to  improve  these  predictions.  Also  note  the  terms  that 
propagate  diagonally  across  the  plate  have  not  been  implemented  at  this  time  and, 
therefore,  are  not  present  in  the  calculated  data. 
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Figure  5.1:  The  2  ft.  square  plate,  tilted  by  45°  with  respect  to  the  horizontal  plane, 
used  in  the  time  domain  analysis;  (a)  Geometry,  (b)  Possible  ray  mechanisms. 
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Figure  5.2:  Impulse  response  of  a  24"  square  plate  at  0  =  30°  and  <p  =  90°.  Vertical 
polarization,  (a)  Measured,  (b)  Calculated. 
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0.003 


Figure  5.3:  Impulse  response  of  a  24"  squa’e  plate  at  ^  =  30°  and  4)  =  90°.  Horizontal 
polarization,  (a)  Measured,  (b)  Calculated. 
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Figure  5.6:  Frequency  response  of  the  edge  wave  diffraction  mechanism  E2  of  a  24" 
square  plate  B.t  6  =  30°  and  <(>  =  90°.  (a)  Vertical  Polarization,  (b)  Horizonte! 
Polarization.  Solid  line:  measurements,  dotted  line:  calculations. 
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Figure  5.8:  Frequency  response  of  the  diffraction  mechanisms  within  the  time  win 
dow  (2ns, 3ns)  of  a  24"  square  plate  at  5  =  75°  and  =  90°.  (a)  Vertical  Polarization 
(b)  Horizontal  Polarization.  Solid  line:  measurements,  dotted  line:  calculations. 


Chapter  6 

Ultra-Wideband  Imaging  of 
Scattering  Centers 

6.1  Goals 

One  of  the  goals  of  the  project  was  to  perform  fundamental  diagnostic  scattering 
mechanism  research  utilizing  both  analysis  and  measurements.  This  goal  was  fulfilled 
using  the  technique  of  ultra- wideband  polarimetric  imaging. 

6.2  Method 

The  technique  of  ultra- wideband  radar  imaging  is  described  in  Report  721198-2,  en¬ 
titled  “Ultra- Wideband  Radar  Imaging  and  the  Diagnosis  of  Scattering  Centers,”  by 
S.  M.  Scarborough.  The  process  uses  measured  or  calculated  complex  backscattering 
data  for  multiple  polarizations  over  a  10:1  range  of  frequency  and  at  least  90  degrees 
of  look  angle  variation.  A  sketch  of  the  “data  space”  is  shown  in  Figure  6.1.  It  is 
known  that  imbedded  in  such  data  is  information  on  not  just  the  location,  but  also 
the  frequency,  angle,  and  polarimetric  “signatures”  of  the  scattering  renters  which 
comprise  the  total  response  of  an  object. 

The  imaging  algorithm  used  in  the  report  is  the  back-projection  technique  [15]. 
In  this  algorithm,  the  image  construction  is  viewed  as  a  weighted  sum  of  one¬ 
dimensional  projections,  each  of  which  is  obtained  from  the  complex  spectral  re- 
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ULTRA  WIDEBAND  RADAR  TYPICAL  CONVENTIONAL 

IMAGE  DATA  WINDOW  IMAGE  DATA  WINDOW 


Figure  6.1:  Measurement  data  space  for  UWB  images 


Figure  6.2:  Image  of  bandlimited  impulse  —  1  back  projection, 

spouse  at  a  particular  target  look  angle  lying  in  the  radar  measurement  plane  of 
angles.  Figure  6.2  shows  an  image  for  a  band-limited  impulse  at  just  one  look  angle. 
Figure  6.3  shows  the  image  for  just  two  look  angles  at  90  degrees  to  one  another. 
Finally,  Figure  6.4  shows  an  image  using  the  band-limited  impulse  over  a  180  degree 
range.  This  figure  represents  the  typical  point  spread  function  of  this  imaging  pro¬ 
cess.  Using  frequencies  from  2-18  GHz  and  angle  data  ranges  of  at  least  90  degrees 
or  typically  180  degrees,  0.1  nsec  or  1.5  cm  resolution  is  obtained. 

A  sampling  theorem  exists  which  governs  how  often  in  angle  and  frequency  the 
response  must  be  recorded.  If  D  is  the  maximum  dimension  of  the  target  and  c  is 
the  speed  of  light,  then  the  sampling  rate  in  frequency  must  be 

-  2(1  +  k)D' 

The  factor  k  in  this  equation  is  a  safety  factor,  whose  value  is  usually  0.5  to  1.  Angle 
sampling  depends  on  both  target  maximum  extent  D  and  the  highest  frequency  of 


67 


Figure  6.3;  Image  of  bandlimited  impulse  —  2  back  projections. 


the  interrogating  radar  signal,  /;/ .  The  frequency  sampling  rate  must  be 


Images  have  been  produced  using  measured  data  of  both  simple  and  complex  tar¬ 
gets.  For  example,  Figures  6.5  and  6.6  show  the  images  for  a  tilted  square  conducting 
plate  at  horizontal  and  vertical  polarization,  respectively.  The  image  simulates  a  top 
view  of  the  object  (perpendicular  to  the  plane  of  the  measurements).  The  outline 
of  the  plate  is  clear  in  both  images.  It  looks  rectangular  instead  of  square  because 
the  plate  is  tilted  45  degrees  with  respect  to  horizontal.  However,  in  addition  to  the 
rectangle  outline  are  the  returns  from  several  other  scattering  mechanisms.  The  cir¬ 
cle  in  the  center  is  the  return  from  the  foam  target  support  cylinder.  The  phantom 
edge  in  the  vertical  polarization  image  is  a  multiple  scattering  mechanism  between 
the  front  and  back  edge  of  the  plate.  The  arcs  in  the  images  are  multiple  scat¬ 


tering  mechanisms  created  by  “edge  waves”  and  corner  interactions  on  the  target. 
The  edge  wave  mechanisms,  in  particular,  have  been  very  difficult  to  identify  before 


Figure  6.4:  Image  of  bandlimited  impulse  over  180°  angular  bandwidth 


this  process  (and  these  data  over  such  a  wide  coherent  frequency  and  angle  range) 
became  available. 

Figure  6.7  shows  the  image  of  a  model  MIG  21  using  horizontal  polarization. 
Figure  6.8  shows  an  image  of  an  M35  army  truck,  with  an  outline  of  the  truck 
superimposed.  Both  of  these  images  show  the  outlines  of  the  structure  are  preserved, 
but  higher  order  scattering  mechanisms  are  also  visible.  It  could  be  that  these  higher 
order  scattering  terms  are  more  distinctive  for  identification  than  the  outlines. 

The  results  also  show  there  is  information  in  these  images  in  addition  to  the  geo¬ 
metrical  distribution  of  the  scattering  centers.  If  a  gate  is  placed  around  a  particular 
scattering  center,  then  a  transform  process  may  be  used  to  extract  the  frequency  and 
angle  response  of  just  that  particular  mechanism.  For  example,  Figure  6.9  shows 
the  square  plate  image,  with  the  outlines  of  three  square  gates  which  were  used  to 
separate  individual  mechanisms.  A  frequency-angle  space  image  of  the  response  of 
the  diagonal  corner  interaction  mechanism  of  window  three  is  shown  in  Figure  6.10. 
A  spectrum  for  the  30  degree  look  angle  extracted  from  the  image  in  Figure  6.10 
is  shown  in  Figure  6.11.  The  corner  interaction  demonstrates  the  behavior  as 
predicted  by  theory. 

6.3  Conclusions 

The  conclusions  from  this  report  are: 

•  The  high  resolution  of  the  images  permits  scattering  mechanisms  to  be  resolved 
distinctly. 

•  The  high  frequency  bandwidth  and  angle  range  shows  that  the  scattering  mech¬ 
anisms  have  distinctive  behavior  vs.  frequency  and  angle,  allowing  them  to  be 
identified  as  well  as  located. 
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Figure  6.5:  Image  of  tilted  square  plate,  -5°(0.5°)95°  HP. 
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Figure  6.6:  Image  of  tilted  square  plate,  — 5°(0.5°)95°  VP. 
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Figure  6.7:  Image  of  model  MIG-21  aircraft,  0°(1°)360°  HP. 
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Figure  6.10:  Frequency /angle  response  of  tip-to-tip  diffraction,  HP. 
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Figure  6.11:  Frequency  response  of  tip-to-tip  diffraction  at  20'’  azimuth.  HP. 
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•  Polarimetric  data  show  it  is  the  individual  scattering  mechanisms  which  have  a 
polarization  response  which  is  simple  enough  to  be  predicted  beyond  the  range 
of  measurements,  or  to  be  interpolated  between  measurements. 

•  When  many  of  these  polarimetric  scattering  mechanisms  are  added,  and  possi¬ 
bly  interact  with  each  other,  the  composite  behavior  of  the  target  can  become 
quite  complex. 
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Chapter  7 


Calibration  and  Performance  of 
Radar  Signature  Measurement 
Systems 

7.1  Goals 

Included  in  the  work  statement  cr»  this  project  were  measurements  of  simple  targets 
such  as  corner  reflectors,  dihedrals,  and  other  simple  reflectors.  Such  data  not  only 
give  information  on  simple  scattering  mechanisms,  but  also  are  the  basis  for  using 
such  targets  as  secondary  calibration  standards  for  polarimetric  scattering  measure¬ 
ments.  These  measurements  have  been  performed  in  conjunction  with  a  study  on 
Ultrr.-Wide  Band  scattering  range  calibration,  and  a  separate  technical  report  has 
been  written  on  this  subject. 

7.2  Method 

The  report  which  contains  this  work  is  entitled  “An  Assessment  of  Three  Experimen¬ 
tal  Scattering  Measurement  Systems  Operating  Below  2  GHz,”  by  John  S.  Gwynne, 
Report  721198-3.  The  report  discusses  the  types  of  measurement  errors  encountered 
in  three  separate  RCS  systems,  and  then  describes  calibration  and  range  evaluation 
techniques  appropriate  for  all  cases.  The  report  concludes  with  some  illustrative 
measurements  of  several  targets  on  the  three  systems. 
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The  three  ranges  studied  were: 


•  The  Ohio  State  University  “Big  Ear”  Radio  Telescope,  used  as  a  compact  radar 
scattering  range  (OSURT) 

•  The  Parallel  Plate  Transmission  Line  scattering  range  at  the  Ohio  State  Uni¬ 
versity  ElectroScience  Laboratory  (PPTL) 

•  The  indoor  Compact  Range  at  the  Ohio  State  University  ElectroScience  Lab¬ 
oratory  (OSUCR) 

Range  construction  and  typical  range  instrumentation  for  all  three  of  these  sys¬ 
tems  is  discussed  in  the  report.  Also  described  are  the  most  common  error  sources 
for  each  system. 

The  report  focuses  on  two  performance  functions  which  characterize  these  ranges. 
The  System  Transfer  Function  (SRT)  defines  the  relative  amplitude  and  phase  vs. 
frequency  and  position  in  the  test  zone  for  an  ideal  point  scatterer.  The  Calibration 
Transfer  Function  (CTF)  is  the  inverse  of  the  SRT,  and  represents  the  complex 
function  which  when  multiplied  by  the  raw  data  wiD  give  calibrated  scattering  data. 

The  report  contains  measured  SRT  data  for  all  three  systems.  The  technique 
for  measurement  is  described,  and  the  noise  and  error  sources  are  described  and 
demonstrated.  The  calibration  flow  diagram  shown  in  Figure  7.1  summarizes  the 
process.  The  results  of  Figure  7.2  show  the  benefit  of  including  all  of  the  steps  in 
the  calibration  process  of  Figure  7.1.  The  selection  of  windowing  functions  for  use 
in  this  process  is  also  discussed  in  detail. 

Beyond  optimizing  the  process  for  system  calibration,  the  calibration  target  has 
an  effect  on  calibration  accuracy  on  a  real  system  with  noise  and  clutter.  This  report 
discusses  the  use  of  spheres,  hemispheres,  and  dihedrals  for  calibration,  and  compares 
their  performance  features.  For  example.  Figure  7.3  shows  a  comparison  of  the  CTF 
measured  using  a  6-inch  dihedral  vs.  one  obtained  from  a  12-cm-diameter  hemisphere 
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Figure  7.1:  Flow  diagram  for  the  calibration  of  data  that  has  a  low  signal  to 
noise/clutter  ratio. 
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Figure  7.2;  Example  of  the  benefits  of  the  additional  calibration  steps  in  the  time 
domain. 
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measurement.  As  seen  in  the  figure,  nulls  in  the  calibration  target  response  create 
spike  artifacts  in  the  CTF  frequency  response. 

Finally,  the  report  presents  performance  characteristics  and  target  measurement 
examples  for  the  three  ranges.  For  example,  Figure  7.4  shows  a  waterfall  image 
^  diagram  of  cylinder  field  probing  echoes  for  the  OSURT.  Error  mechanisms  and 

their  relative  amplitude  compared  to  the  desired  plane  wave  beam  are  easily  visible 
in  this  image. 

Figure  7.5  shows  the  outline  diagram  of  an  inert  land  mine  model  which  was 
measured  on  all  three  ranges.  Figure  7.6  shows  its  calibrated  backscattered  echo 
spectrum  at  broadside,  measured  on  the  three  ranges.  Figure  7.7  shows  its  ap¬ 
proximate  impulse  response  derived  from  the  measured  spectral  data.  The  report 
describes  the  factors  which  cause  the  slight  differences  in  these  data. 

7.3  Conclusions 

The  conclusions  section  of  the  report  describes  the  relative  strengths  and  weak¬ 
nesses  of  the  three  measurement  systems.  It  stresses  the  importance  of  choosing  a 
calibration  target  which  has  no  pattern  or  frequency  response  nulls,  and  also  has  an 
accepted  accurate  theoretical  solution  for  its  response  signature.  The  section  also 
discusses  improvements  which  could  be  made  to  the  ranges  and  to  their  calibration 
targets  and  process. 
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Figure  7.4:  Waterfall  diagram  of  cylinder  field  probing  results  for  the  OSURT 
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Figure  7.6:  Measurement  example:  6"  mine  at  broadside-frequency  domain. 
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Figure  7.7:  Measurement  example:  6"  mine  at  broadside-time  domain. 
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Chapter  8 
Conclusions 


This  research  effort  has  sought  to  investigate  electromagnetic  scattering  phenomenol¬ 
ogy  over  a  wide  range  of  frequencies  and  polarizations  to  enhance  performance  of 
SAR  ground  imaging  radars  for  automatic  target  detection/classification.  The  effort 
has  involved  closely  coordinated  theoreticsd  and  experimental  studies,  and  included 
expanded  data  base  measurements  not  envisioned  in  the  original  contract  work  state¬ 
ment. 

Analytical  results  include  the  following: 

•  Theoretical  models  for  dihedral  scattering  signatures  have  been  completed. 

•  Extensions  of  the  dihedral  models  have  been  accomplished  for  the  top-hat, 
bi-cone,  “cake  pan,”  bruderhedral,  and  “donut”  shapes. 

•  An  improved  far  zone  diffraction  coefficient  for  a  vertex  of  an  angular  sector 
has  been  developed. 

•  Edge  wave  and  double  diffraction  mechanisms  across  a  plate  have  been  studied 
to  determine  their  frequency  and  pattern  behavior  in  images. 

•  Preliminary  double  diffraction  coefficients  for  non-coplanar  skewed  edges  nec¬ 
essary  for  the  trihedral  analysis  have  been  developed. 

Experimental  results  include  the  following: 
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•  UWB  polarimetric  scattering  signatures  of  real  and  scale  model  targets  were 
measured. 

•  UWB  diagnostic  backscattering  images  were  generated  from  measured  data 
for  both  canonic  and  realistic  targets. 

<  i 

•  A  UWB  calibration  procedure  was  developed,  and  was  developed  on  three 

different  systems.  ^ 

•  A  data  processing  technique  to  extract  the  frequency  spectrum  and  pattern 
behavior  of  individual  scattering  centers  from  UWB  diagnostic  images  was 
developed  and  demonstrated. 

•  The  speed,  dynamic  range,  and  accuracy  of  the  “Big  Ear”  compact  range 
radar  system  were  improved.  The  system  was  useful  for  targets  as  large  as  70 
ft.  long. 

In  conclusion,  the  studies  have  shown  that  there  do  exist  scattering  mechanisms 
which  are  strongly  related  to  target  geometry  and  hence  are  robust  over  polarization, 
frequency,  and  look  angle.  We  have  established  the  procedures  for  determining  a 
wide  variety  of  such  scattering  centers  for  a  realistic  object,  either  analytically  or 
experimentally.  We  recommend  that  these  scattering  centers  be  prime  candidates 
for  target  discrimination  research  in  the  future.  We  also  recommend  further  studies 
of  the  behavior  of  such  mechanisms  for  targets  in  foliage  and  other  complex  clutter- 
limited  situations. 
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